Classical and quantum interference in multiband optical Bloch oscillations 
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Classical and quantum interference of light propagating in arrays of coupled waveguides and un- 
dergoing multiband optical Bloch oscillations (BOs) with negligible Zener tunneling is theoretically 
investigated. In particular, it is shown that Mach-Zehnder-like interference effects spontaneously 
arise in multiband BOs owing to beam splitting and subsequent beam recombination occurring in 
one BO cycle. As a noteworthy example of quantum interference, we discuss the doubling of inter- 
ference fringes in photon counting rates for a correlated photon pair undergoing two-band BOs, a 
phenomenon analogous to the manifestation of the de Broglie wavelength of an entangled biphoton 
state observed in quantum Mach-Zehnder interferometry. 

PACS numbers: 72.10.Bg, 72.90.+y, 42.82.Et, 42.50.Dv 



I. INTRODUCTION 



Since the pioneering works by Bloch, Zener and 
Wannier on the coherent electron dynamics in biased 
crystalline solids [H, it is well established that quantum 
particles in periodic potentials subjected to an external 
force do not delocalizc but undergo a high-frequency 
oscillatory motion known as a Bloch oscillation (BO). 
Tunneling to higher-order bands, referred to as Zener 
tunneling (ZT) and observed for a strong external 
bias, is responsible for BO damping and broadening 
of Wannier-Stark resonances. Owing to detrimental 
dephasing and many-body effects, the experimental 
verification of BOs failed for many decades. Nowadays 
BOs are considered a rather universal wavy phenomenon 
and their observation has been reported in different 
physical systems for both quantum particles and classical 
waves. In the quantum realm, BOs have been observed 
for electrons in biased semiconductor superlattices 
and for Bose-Einstein condensates in accelerating 
optical lattices [2>\. Classical analogues of BOs, based on 
interference effects of classical waves in periodic media, 
have been proposed and experimentally observed using 
either optical 0-0| or acoustical [|| waves. Optical BOs 
in artificial materials or at the nanoscale have been 
also recently predicted @. Among various photonic 
structures, arrays of coupled waveguides with transverse 
refractive index gradients have provided an extremely 
rich laboratory tool to visualize the classical wave optic 
analogues of BOs [H H, 0] and related phenomena, such 
as the coherent superposition of BOs and ZT (the so- 
called Bloch-Zener oscillations [T3|) occurring in binary 
structures and recently observed in circularly-curved 
femtosecond-laser-written waveguide arrays 111. In 
another physical context, application of Bloch-Zener 
oscillations to matter-wave interferometry has been also 
proposed [12]. While there has been a lot of interest 
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in the propagation of light waves in complex photonic 
structures like lattices and superlattices, quasi-crystals, 
and even random lattices, most of previous works focused 
on propagation of classical light and used the photonic 
structures as 'classical simulators' of electronic or matter 
wave systems in periodic or random potentials. However, 
it is known that coupled waveguides behave similarly to 
beam splitters (see, for instance, [1J]) and may therefore 
show quantum interference effects when probed with 
nonclassical light. Demonstration of key quantum effects 
of nonclassical light in silicon-based based waveguide 
circuits has been recently reported in Ref . [lj] . As shown 
in recent works [HI, [lg|, nonclassical light consisting 
of only particle-like quanta propagating in waveguide 
arrays with a superimposed transverse refractive index 
gradient can also produce optical BOs. Remarkably, in 
addition to classical wave Bragg scattering the quantum 
nature of light introduces new quantum interference 
effects. In particular, two-photon Hong-Ou-Mandel 
quantum interference 17[ has been predicted for pairs of 
correlated photons undergoing Bloch-Zener oscillations 
in binary arrays. Here ZT periodically mixes photons 
belonging to the two minibands of the array and thus 
acts like a beam splitter [l5j . Recently, quantum corre- 
lations of photon pairs undergoing discrete diffraction 
in homogeneous arrays have been investigated as well, 
and their classical counterpart has been experimentally 
observed in Hanbury Brown-Twiss intensity correlation 
measurements (l8j . The investigation of quantum 
interference effects in complex photonic structures 
such as homogeneous, inhomogeneous or even random 
photonic lattices has two main motivations. On the 
one hand, they may offer the possibility of engineering 
photon entanglement and of transporting nonclassical 
light [l8|, [l9| ; on the other hand, they enable to utilize 
optical-quantum analogies in the opposite direction, 
i.e. trying to observe quantum version of classical wave 
phenomena (such as the analogue of the Talbot effect in 
the second quantized setup |19j). 

It is the aim of this work to provide a comprehensive 
analysis of classical and quantum interference of light 
undergoing BOs in inhomogeneous waveguide arrays 
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in the regime of multiband excitation and negligible 
ZT, i.e. in a complementary regime of that previously 
considered in Refs.JH, HI]- In spite of the absence of 
interband transitions, it is shown that Mach-Zehnder- 
like interference effects do spontaneously occur for a 
multiband excitation of the array at the input plane. 
In this regime, the input beam breaks up into two 
(or more) wave packets belonging to different bands 
of the array which follow distinct paths in the real 
space and recombine after a full BO cycle. Wave and 
photon interference effects after a full BO cycle may 
be thus observed by classical and nonclassical light 
illumination, even in absence of ZT. In particular, it 
is shown that array excitation with correlated photon 
pairs tilted at the Bragg angles enables to observe 
quantum interference patterns in two-photon correlation 
measurements analogous to those occurring in quantum 
Mach-Zehnder interferometry as a manifestation of the 
de Broglie wavelength of entangled photon states |20l - |23l | . 



II. WAVE-OPTICS MODEL OF MULTIBAND 
BLOCH OSCILLATIONS IN A WAVEGUIDE 
ARRAY 

The starting point of our analysis is provided by 
a rather standard wave optics model describing BOs 
of monochromatic light waves at carrier frequency 
to = 27rco/A propagating in a weakling guiding one- 
dimensional waveguide array [Fig. 1(a)], with a periodic 
refractive index profile n{x) and with a superimposed 
transverse refractive index gradient Fx @, [H, H3|- I n 
the paraxial approximation, the slow evolution of a scalar 
field component 4>(x, z) along the paraxial z direction is 
governed by the Schrodinger-like wave equation (see, for 
instance, |2J|) 
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FIG. 1: (color online) (a) Schematic of a singly-periodic one- 
dimensional waveguide array, (b) Behavior of the array refrac- 
tive index V(x) = n 3 — n(x) used in numerical simulations. 
The transverse index ramp term is applied by e.g. transverse 
thermal heating or by circularly-curving the waveguide axis 
(not shown in the figure), (c) Band diagram, (d) Behavior of 
the Fourier coefficients \B n (k — K + 2irl/a)\ 2 = |#n,i(At)| 2 of 
Bloch modes </?„(x, k) for a few low-order bands. 



i*<t>* = 7) $xx + [V(x) - Fx](t>, 
2n s 



(1) 



where X = X/(2tt) is the reduced wavelength, n s is the 
substrate refractive index, V(x) — n s — n(x) accounts for 
the periodic modulation of the refractive index with spa- 
tial periodicity a [V(x + a) = V(x)], and F is the super- 
imposed transverse index gradient. The array is assumed 
to be confined in the region < z < L [see Fig. 1(a)] and 
array excitation is accomplished at the input plane z — 0. 
Since V(x) = F = for z < and z > L, usual paraxial 
diffraction occurs before and after the light beam enters 
and leaves the arrayed structure. The sample length L is 
typically chosen to be equal to the characteristic period 
zb of BOs at a reference value F = F of index gradient, 
i.e. L — zb where [23 



zb 



Fa' 



(2) 



In absence of the transverse index gradient (F = 
0), the modes of the arrayed structure are the Bloch 



states ip n (x,K,) with corresponding dispersion curves 
E n (n), where n = 1,2,3,... is the band index, n 
is the Bloch wave number (quasi momentum) which 
varies in the first Brillouin zone —ir/a < k < 
7r/a, and the orthogonal and normalization conditions 

(tfn' {.X, K,')\ip n {x, K)) = 5 n}n '5(K - «') hold. (f n (x, k) 

and E n (n) are found as eigenfunctions and eigenvalues 
of the problem 'Hof n (x 1 k) — E n («)</?„ (x, k) with pe- 



riodic Hamiltonian Hq 



-[X'/(2n s )]d 2 xx + V(a 



Ac- 



cording to the Bloch theorem, one can write <p n (x, k) = 
u n (x, k) exp(ira), where the periodic part u n (x, n) of the 
Bloch state can be expanded as a Fourier series 



u n (x, k) 



1 



DO 



2tt , 



l ,i(/v)exp[(2i7rZ/a)x] (3) 



with coefficients n ,i{^)- Owing to the orthogonal and 
normalization conditions of Bloch states, the matrix 
A = 0n,i turns out to be unitary, i.e. A^ 1 = A'. In 
place of the Fourier coefficients 6 n .i(K), defined in the 
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first Brillouin zone — n/a < k < n/a, one can introduce 
the set of functions 



B n {k 



2nl/a) 



(4) 



which are defined over — oo < k < oo. A typical ex- 
ample of refractive index profile and corresponding band 
diagram is shown in Figs. 1(b) and 1(c) for an excitation 
wavelength A = 633 nm, lattice period a = 11 fim, bulk 
refractive index n s = 1.43 and maximum index change 
of the waveguides An = 0.0016. The behavior of the 
Fourier coefficients |i? n (fc)| 2 of Bloch states is also de- 
picted in Fig. 1(d). 

In presence of a transverse index gradient F, small 
enough to neglect inter band transitions (i.e. ZT), the so- 
lution to Eq.(l) for an assigned beam distribution 4>{x, 0) 
at the z — input plane can be written as a superposi- 
tion of wave packets <fi n (x,z) belonging to the different 
bands of the arrays which propagate independently each 
other and undergo BOs, namely one can write [24| 



(x,z)= ^2 ^n{x,z). 



(5) 



n=l,2,3, 



For negligible ZT, the evolution of the wave packet </>„ can 
be calculated in a closed form by using the acceleration 
theorem [25[ and reads [24| 

<-7r/a 

<f> n (x,z)= / dKCn(K)<p n (x,K+Fz/X)exp[—i'y n (K,z)]. 

J — ir / a 

(6) 

In Eq.(6), the coefficients c n (K,) are determined by the 
field distribution illuminating the array at the input 
plane according to 



c n (/c) = ((p n (x,K)\<j>(x,0)} 



dxipl(x,K)(l)(x,0) (7) 



whereas the phase term j n (n, z) is given by 
I rFz/X 

Jn(K,z) = — dn'E' n (K + «'), 



F 



where 



2tt f a 

E' n (n) = E„(k) - F — / dx iu n (x, 
a Jo 



K) 



du n (x, k) 
dn 



(8) 



(9) 



is the band dispersion curve, corrected to account for a 
possible topological (Berry phase) contribution [26j, l2JJ . 
In particular, after a full BO cycle, i.e. for a propagation 
distance z — zb, one has <j) n {x, zb) = 4>n( x , 0) exp(— vy n )i 
where 



1 



2-ir/a 



dK,E n {n). 



(10) 



Since in the general case the phases 7„ accumulated by 
the various wave packets <f> n undergoing BOs are not the 
same or do not differ by multiplies of 2tt, the output field 



distribution 4){x,zb) does not generally reproduce the 
input one <j)(x,0). The output field is in turn given by 
the interference of the various wave packets (f> n (x, 0) with 
the appropriate phase delays j n given by Eq.(10). This 
kind of interference is discussed in details in the next sec- 
tion and may lead to a Mach-Zehnder-like interferometry 
which does not require ZT. 



III. BEAM BREAK UP AND 
RECOMBINATION: CLASSICAL 
INTERFERENCE 

Let us consider a broad input beam, with narrow an- 
gular spectrum and near-field distribution G(x) , imping- 
ing the array at the incidence angle 9, so that (f>(x, 0) = 
G(x) exp(iKx) where k = n s 9/X = {jt / a){6 / 6 b) and 



A 



IB = 



2cm, 



(11) 



is the Bragg angle [281] . We denote such an input field 
distribution as (x, 0; Ko) and the corresponding propa- 
gated field as g^(x, z; no), where kq and the integer I are 
defined such that k + 2irl/a = k and —n/a < k < ir/a. 
The normalization condition J dx\g^ (x, z; Ko)\ 2 = 1 is 
assumed. In this case, following the analysis of Ref. [HJ , 
one can show that inside the arrayed structure one has 



\0n(x,z)\ 2 

where we have set 



n,j(Ko)¥»n \ X, Kq + ^-^J G(x - X n (z)) 

(12) 



x n {z) = 1 [E' n ( Ko + Fz/X) - E'(k )} . (13) 

Equations (5) and (12) indicate that the injected beam 
4>(x, 0) breaks into a superposition of wave packets <f> n be- 
longing to the different bands of the array, with weight- 
ing factors \0 n j(Ko)\ 2 , which undergo BOs along different 
paths x n (z) [see, for instance, Figs. 3(a) and 4(a)-(d) to be 
discussed below]. According to Eq.(13), the path x n (z) 
followed by the wave packet tp n reproduces the shape of 
the band dispersion curve E n (k) , eventually corrected to 
include the Berry phase contribution. After a full BO 
cycle, i.e. at z — zb, the different wave packets </>„ in- 
terference and the following scattering relations hold (see 
the Appendix) 

oo 

g {l) (x,z B ;Ko) = Mi, P g ip) (x,0; Ko ) (14) 

p— — oo 

where the scattering matrix A4 is given by 

M = A^BA, A n ,i = 0»u(«o), B P ,n = exp(-i'yn)&p,n, 

(15) 

and 7„ are given by Eq.(10). Note that, if the differ- 
ences of phase delays 7„ — 7 P were multiplies of 2-7T, a full 
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reconstruction of the input beam -apart from an unim- 
portant phase term- would be achieved since A^ A = I. 
In this case the beam leaving the array at z = zb would 
not split and would diffract at the same tilting angle 9 as 
that of the incoming beam. However, the phase delays 
7„ do not generally satisfy the previous condition, and 
thus the beam leaving the array at z — zb breaks up 
into several beams which diffract at the different angles 
9 + 2p6 B (p = 0,±1,±2,...) according to Eq.(14). Af- 
ter some propagation distance, such beams are not over- 
lapped and can be thus spatially resolved. This general 
behavior is schematically illustrated in Fig. 2(a). More 
generally, we can say that, if the array is simultaneously 
excited by a set of identical broad beams tilted at differ- 
ent angles 9 + 210b (I = 0. ±1. ±2. ...), the array behaves 
like a liner multiport optical system for the amplitudes 
of modes [see Fig. 2(b)] with a scattering matrix M. 
given by Eqs.(15). Here we have focused our analysis 
to the case where the length L of the array equals the 
BO cycle zb- Similar results are obtained by assuming 
-more generally- that the array length L comprises an 
integer number N of BO cycles, i.e. for L = Nzb- In 
this case, the scattering matrix can be again factorized 
as M. = A^BA, where now B p . n — exp(—iNj n )6 p>n . 
As an important example, let us consider the case where 
only the two lowest bands of the array are involved in 
the dynamics, and the scattering relations (15) may be 
limited to two modes solely. An inspection of the curves 
E n {n) and \B n (k)\ 2 shown in Figs.l (c) and (d) indicates 
that this condition is realized, to a good approximation, 
when the array is excited by two broad beams tilted the 
former at an angle 9\ equal to or slightly smaller than the 
Bragg angle 9b, the latter at the angle 92 = 9\ — 29b ~ 
—9b [see Fig. 3(a)]. If we use the simplified notations 
gi(x, z) = gW(x,z;K ) and g 2 (x,z) = z; k ) for 

the two modes scattered by the array in one BO cycle, 
where kq ~ 7r/a, the input-output relations [Eq.(14)] re- 
duce to 

gi(x,z B ) = Sngi(x,0) + Si 2 g 2 (x,0) (16) 
g 2 (x,z B ) = S 2 igi(x, 0) + S 22 g 2 (x, 0) (17) 

where we have set Su = M.0,0, Su = A^o,-i, S 2 i = 
M-ifi and S 22 = M.-1-1. According to Eqs.(15), the 
2x2 scattering matrix S entering in Eqs.(16) and (17) 
is given by 

s _ ( P11 P21 \ x ( exp(-ni) \ x / p u p u 
V Pu Ph J V exp(-i7 2 ) J \ P21 p 22 

(18) 

where 

( Pu Pu \ = ( 61,0(1*0) °i-i( K o) \ , ig x 

V P21 P22 J \&2,o(Ko) ^2,-l(«o) /' 

Since the matrix pu- is unitary, one has \pn\ = | P22 1 > 
\pu\ = |p2i|, |/On| 2 + |pi 2 | 2 = 1 and pup* 2 + p 2 ip 22 = 0. 



Without loss of generality, we may assume pu = p 22 to 
be real-valued and positive by a suitable choice of the 

(a) (b) 




x 



FIG. 2: (color online) (a) Schematic of BO motion in the 
multiband regime. A broad beam illuminates the array at 
the incidence angle 9 and excites several wave packets 4>i , 4>i, 
<j>3 belonging to the different bands of the array. The wave 
packets undergo BOs along distinct paths which reproduce in 
real space the spectral shapes of the band dispersion curves. 
After one BO cycle, owing to different phase delays 7„ accu- 
mulated by the wave packets, several beams, diffracted at the 
angles 6, 6 ± 26b, 9 ± 4#b, are produced, (b) Schematic 
of a linear-optic multiport system described by a scattering 
matrix M . 



absolute phases of u\(x, kq) and u 2 (x, kq). After setting 
pu = P22 = VT and R=l — T, where T = |6»i, M| 2 = 
\@2.-i(ko)\ 2 , we can then write 

f Pu Pu \ _ ( _VT v A Rexp(ia) \ 

V P21 P22 ) V — /Rexp(-ia) VT J 

where a is the phase of pi 2 . Note that Eq.(20) is analo- 
gous to the scattering matrix of a lossless beam split- 
ter with transmittance T = 1 — R (see, for instance, 
[29} ). Physically, such a transformation corresponds to 
the mixing of the incoming beams g\ and g 2 into the 
two wave packets 4>i and 4> 2 belonging to the two lowest- 
order bands of the array [see Fig. 3(a)]. At exact Bragg 
incidence, i.e. for 6\ = —9 2 = 9b, one has T = 1/2, i.e. 
the equivalent beam splitter is balanced. Note also that 
the full transformation (18) is analogous to that of a two- 
port Mach-Zehnder interferometer in which two waves g\ 
and g 2 are mixed by a fist beam splitter BS1, undergo 
different phase delays 71 and 72 in the two arms of the in- 
terferometer, and are then recombined by a second beam 
splitter BS2 [see Fig. 3(b)]. The transfer matrix S of the 
equivalent Mach-Zehnder interferometer reads explicitly 
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( Sn S12 \ _ ( Texp(-Z7i) + i?cxp(-z7 2 ) \/RTexp(ia) [exp(-i7i) - exp(-i7 2 )] 

\ S21 S22 J \ \J RT exp(— ia) [exp(— 171) — exp(— 1^2)] Pexp(— 272) + Pexp(— 271) 



Classical Mach-Zehnder interferometry, which do not re- 
quire ZT, is thus expected to be observable in multiband 
BOs. As an example, for single beam excitation at the 
tilt angle 0\, i.e. for (f>(x, 0) = g\ (x, 0), the fractional light 
powers Pi(A7) and P 2 (A7) of the two beams gi(x, zb) 
and 52(^1 z b) leaving the array and measured by two pho- 
todetectors Dl and D2 [see Fig. 3(a)], are readily calcu- 
lated from Eqs.(16), (17) and (21) as 

Pi(A 7 ) = l^nl 2 = T 2 + i? 2 + 2i?Tcos(A7) (22) 
P 2 (A 7 ) = |S 2 i| 2 = 2PT(l-cos(A 7 )). (23) 

According to Eq.(10), the phase difference A7 is given 
by 

2tt / a [ 27T/a , r, M , . .t\ 2irAE 
where AP is the distance between the first and scc- 




x 



FIG. 3: (color online) (a) Schematic of two-band BOs in a 
waveguide array excited by two input broad beams tilted at 
the angles 81 ~ 6b and 62 — 0i — 29b ~ —Ob- The paths 
followed by the wave packets <j>\ and 02 undergoing BOs map 
in real space the spectral shapes of the band dispersion curves 
Ei(k) and Ez(k), respectively, (b) Equivalent two-port Mach- 
Zehnder interferometer. 

ond bands of the array, measured from their dc values 
(27| . Typically the value AP is of the order of the refrac- 
tive index change An of the core region of waveguides 
from the dielectric substrate [see Fig. 1(c)]. In practice, 
to vary A7 one can slightly change the index gradient 
P around the reference value Pq = X/(aL) [3(j. In this 



I 

way, the deviation of Zb = A/ (Pa) from L = \/(Fqo) 
is negligible, i.e. the BO cycle is almost completed at 
the output of the arrayed region, whereas the change of 
A7 can be of the order of 2tt or larger as Fa is typi- 
cally much larger than AP. Figures 4(a-d) show typical 
examples of two-band BOs observed in numerical simu- 
lations of Eq.(l) for the one-dimensional array of Fig.l. 
The length of the arrayed region is L — 23 mm, and the 
condition zb = L is attained for an applied refractive 
index gradient Pq = 2.502 m" 1 . For such a relatively 
low value of refractive index gradient, ZT from band 1 to 
band 2, and from band 2 to band 3, turns out to be neg- 
ligible. The array is excited by a broad Gaussian beam 
tilted at the Bragg angle 62 = —&b, and the evolution 
of beam intensity \4>(x, z)\ 2 along the sample is plotted 
for a few values of the refractive index gradient P close 
to Fq. Beam splitting and beam recombination after a 
full BO cycle are clearly visible, as well as the change of 
the power levels in the two scattered output beams as 
the applied index gradient P is slightly varied (by a few 
percents) at around Fq. The behavior of the fractional 
powers Pi and P2 carried by the two output beams versus 
the ratio F /F is depicted in Fig. 4(e), clearly showing an 
oscillatory behavior with a visibility of about ~ 90%. 

IV. QUANTUM INTERFERENCE IN 
MULTIBAND BLOCH OSCILLATIONS WITH 
NONCLASSICAL LIGHT 

In previous sections, multiband optical BOs and in- 
terference effects have been investigated in the frame- 
work of classical electromagnetic theory. Here we ex- 
tend the previous analysis to nonclassical states of light 
undergoing multiband BOs, and discuss relevant quan- 
tum interference effects similar to those observed in 
Mach-Zehnder-based quantum interferometry [2"oT - [23l | . A 
second-quantization approach to study optical BOs in 
waveguide arrays, based on a procedure similar to that 
adopted in the quantum theory of solitons [3l| , has been 
developed in Ref.[l5|]. In this approach, the scalar wave 
equation for the classical field (f>(x,z) [Eq.(l)] is written 
in Hamiltonian form assuming the paraxial spatial coor- 
dinate z as an independent variable, and a quantization 
procedure is then applied by introducing creation (j>'(x) 
and annihilation <p(x) bosonic field operators (for details 
see [15I HH). In the Schrodinger picture, the quantum 
field is described by a quantum state \Q(z)) which can 
be expanded in Fock space as 

00 

|Q(z)) =]>>„!/(")), (25) 

n=l 
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FIG. 4: (color online) Classical Mach-Zehnder interferometry 
based on two-band BOs. (a),(b),(c) and (d) show the evolu- 
tion of light intensity distribution along the optical structure, 
excited at the input plane by a broad Gaussian beam tilted 
at 62 = —0b, for a few values of Fo/F. The arrayed structure 
has a length L — 23 mm and it is comprised between the 
two dotted horizontal lines. The index gradient Fo, at which 
the array length L is exactly equal to one BO cycle zb, is 
given by Fo = 2.502 m _1 . In (e) the numerically-computed 
fractional powers Pi and P2 carried by the two output beams 
(right- and left-diffracted beams, respectively) are plotted as 
functions of the ratio Fo / F. 



where the n-photon number state |/^"^) is defined by 



|/(«)) = / d Xl dx 2 ..Ax / n) ^ X Jl-" Xn ' z) 
x <pt(x 1 )^(x 2 )...^(x n )\0). 



(26) 



The normalization conditions \a n \ 2 = 1 and 
/ dx\dx2---dx n \f^ (xi, a?2, x n , z)\ 2 = 1 are assumed. 
The Fock state |/^) defined by Eq.(26) is obtained from 
the vacuum state |0) by creating n photons with a spa- 
tial weighting function /(™)(xi, X2, x n , z). The evo- 
lution equation for obtained from the Schrodinger 
equation and us ing the commutation relations of field 
operators, reads [151 ] 



iX 



df(n 

dz 



E 

1=1 



-kh +v{xl) ~ Fxi 



f 



in) 



(27) 



Owing to the bosonic nature of photons, solely symmet- 
ric functions should be considered. The simplest n- 
photon number state, denoted by \g) n , is obtained by as- 
suming = g(xi, z)g(x2, z)...g(x n , z), where the func- 
tion g(x, z) satisfies the classical wave equation (1) with 
the normalization J dx\g(x, z)\ 2 = 1. Note that in this 
case one has 



\g) n = -;=}[ dxg(x,z)ip(x)) |0), 



(28) 



so that this quantum state describes the excitation of 
the array with a beam with a spatial profile g(x,0) at 
the entrance plane z = and carrying n photons. More 
generally, for a given set of orthogonal and normalized 
functions gi(x, z), g 2 {x, z), ... that satisfy Eq.(l), one can 
construct the n-photon number state \gi, g2> ■■■)m,n 2 ,— 
defined by 



"2,- 



dxgi(x, z)<fi^(x) I x 



dxg 2 (x, z)(j)^(x) I x ....|0) 



(29) 



with ni+ri2 + ... = n. Such a quantum state describes the 
excitation of the array with a set of independent beams 
each carrying n\, n.2,... photons (n = m + ri2 + ...). The 
classical picture of multiband BOs, described in previ- 
ous sections, is retrieved from the quantum model when 
the input beam is in a coherent state (classical light), i.e. 
when a superposition of Fock states with a Poisson distri- 
bution is considered [l5| ■ Here we focus to the case where 
the input beam describes a nonclassical field. In partic- 



ular, let use suppose that the array is excited by a set of 
tilted beams g"'(x, 0; Ko) (I = 0, ±1, ±2, ...), which have 
been introduced in Sec. Ill, and assume that is in a 
photon number state (a Fock state) carrying ni photons. 
After introduction of the creation operators a\ 



at = / dxg il) (x,0:,K O )(j>Hx), 



(30) 
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the quantum state of the system at the input plane of 
the array may be written as 



|G(* = o)> = 



1 



- no ~ ni ~"2 



•v/...n !ni!n 2 !. 
|...,n ,ni,n 2 , 



•10} 



(31) 



Note that, as <7^(a;, 0; «o) are orthogonal functions, the 
operators a\ and a/ satisfy the commutation relations 



[at, a 



ti - 



>l,p , l a i,a p \ = [al , a 



tl - 



(32) 



The quantum state of the system at the output plane, 
i.e. after one BO cycle, can be derived from Eqs.(14) 
and (29), and reads explicitly 



|fi(* = *u)) 



1 



^/ ...n \ni\n 2 \ 



r 



M 2 , P a p 



•10). 



(33) 



Note that, since the multiband BOs problem over one 
BO cycle admits of a formulation in terms of a scattering 
matrix [Eq.(14)], the quantum state at the output plane, 
as given by Eq.(33), is consistent with the result that one 
would obtain using the input-output operator formalism 
commonly adopted for linear quantum-optical networks 
(see, for instance, (33[)- Equation (33) is at the basis of 
quantum interference and quantum entanglement observ- 
able in multiband BOs when the array is excited by pho- 
ton number states. As an example, we discuss in detail 
the doubling of interference fringes in photon counting 
rates for correlated photon pairs undergoing two-band 
BOs, a phenomenon analogous to the doubling of inter- 
ference fringes in photon correlation measurements ob- 
served in a Mach-Zehnder interferometer and attributed 
to the so-called de Broglie wavelength of the entangled 
biphoton state (see, for instance, 22]). To this aim, let us 



J 



consider the two-band BO problem of Fig. 3(a) discussed 
in Sec. III. Instead of using classical light (i.e. coherent 
states for the two beams g± and 32), let us illuminate the 
array by correlated photon pairs generated by frequency- 
degenerate spontaneous parametric down-conversion and 
incident onto the array at the Bragg angles Q\ = 9 b and 
82 = —Ob- The quantum state of light at the entrance 
plane of the array can be thus written as 



|(2(z = o)) = |i)i|i> 3l 



(34) 



where |n,i)i|n.2)2 denotes a n = n\ + 712 photon num- 
ber state with n\ photons in the mode g\ and n 2 pho- 
tons in the mode <?2- Using Eqs.(16), (17) and (34), from 
Eqs.(31) and (33) it readily follows that the state of quan- 
tum field after one BO cycle is given by 



\Q(z = z B )) = >/2SnS2i|2)i|0>2 + {S11S22 + S12S21) |l>i|l>2 + >/2Si 2 5 2 2|0) 1 |2)s 



(35) 



r 



The joint probabilities i?n and -R22 to find both photons 
in the same beam, either g\ or g 2 , after one BO cycle are 
then given by 



R%1 — 2|<Sii| 152 



i?2 



2 1 Sis 



(36) 



Such relations give the two-photon counting rates that 
one would measure in an experiment (22j . Using Eq.(21) 
and assuming T = R = 1/2 (valid for incidence angles 
exactly tuned at the Bragg angles ±# B ), one finally ob- 
tains 



R11 = R2 



1 



(1 - cos(2A 7 )) . 



(37) 



Note that the counting rates Rn = R22 oscillate like 
~ cos(2A7), i.e. at twice the phase difference A7, as op- 
posed to the classical first-order interference fringes [3J] , 
which oscillate like ~ cos(A7) [see Eqs.(22) and (23)]. 
The doubling of the counting rate oscillation frequency 
versus the phase delay A7 is analogous to that observed 
in an ordinary Mach-Zehnder interferometer probed by 
correlated photon pairs and generally explained as a man- 
ifestation of the Broglie wavelength of the biphoton en- 
tangled state formed after the first beam splitter BS1 and 
probed at the second beam splitter BS2 of the interfer- 
ometer 



22] 
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V. CONCLUSIONS 

In this work interference phenomena for classical 
and non-classical light propagating in arrays of coupled 
waveguides and undergoing multiband optical Bloch os- 
cillations with negligible Zener tunneling have been the- 
oretically investigated. A wave scattering analysis of 
multiband BOs shows that Mach-Zehnder-like interfer- 
ence effects spontaneously arise owing to beam splitting 
and subsequent beam recombination occurring at each 
BO cycle. A noteworthy example of quantum inter- 
ference is provided by the doubling of the interference 
fringes in photon counting rates for a correlated pho- 
ton pair undergoing two-band BOs. This phenomenon is 
analogous to the one observed in a Mach-Zehnder inter- 
ferometer excited by pairs of correlated photons and is a 
manifestation of the so-called Broglie wavelength of the 
two-photon entangled state produced after the first beam 
splitter of the interferometer and probed by the second 
one [22j]. It is envisaged that the present results may 
stimulate further theoretical and experimental investiga- 
tions of classical and quantum interference phenomena of 
light propagating in complex periodic, quasi-periodic or 
disordered photonic lattices. 

Appendix A: Derivation of the scattering matrix 

In this appendix we derive the scattering relations 
given in the text [Eqs.(14) and (15)] for the amplitudes 
g") of tilted beams at the input and output planes of 
the array, i.e. after one BO cycle. To this aim, let us 
first prove Eqs.(14) and (15) when the incident beams 
are tilted plane waves, i.e. let us first assume 

(x, 0; Ko) = exp ( in^x + i — Ix j (Al) 

V27T V a J 

(I = 0, ±1, ±2, ...). Owing to the completeness of Bloch 
states fn(x, k), one can write 

/it/a 
dn{Lp n \g { - l) )ip n {x, n). 
n=±,*.o,... -^/ a 

(A2) 

Since (f n (x,n) = u n (x, k) exp(iKir) and taking into ac- 
count the Fourier decomposition of u n (x, n) [Eq.(3)], one 



readily finds ((p n (x, n)\g">(x, 0; k)) = 0* n l {n Q )8{K- k ), 
so that Eq.(A2) yields 

g {l) (x,0- 1 K O )= A,nVn{x,K Q ) (A3) 

n=l,2,3,... 

where we have introduced the matrix A n .i = n ,i{i^o) and 
A^ is the adjoint of A, i.e. A] n = A* nl . After one BO 
cycle, the Bloch state ip n (x, Ko) accumulates a phase shift 
exp(— i^ n ), where 7„ is given by Eq.(10). After the in- 
troduction of the diagonal matrix B pM = cxp(— vy n )$n,p, 
we can thus write 

g il) (x,z B ;K Q ) = ^2 Al n cxp(-ij n )if n (x,K ) 

n=l,2,3,... 

= J2 (■4 t B)i,»V»(s.«o)- (A4) 

n=l,2,3,... 

Since the matrix A is unitary, A^ 1 = A' and Eq.(A3) 
can be inverted yielding 

lfi n (x,K Q )= 2J An^g^Hx^;^). (A5) 

;=o,±i,±2,... 

Substitution of Eq.(A5) into Eq.(A4) finally yields 
g^(x,z B ; Ko )= ]T (^BA)i, p g^{x,0;K ) 

p=0,±l,±2,... 

(A6) 

which are precisely the scattering relations between in- 
put and output waves given in the text [Eq.(14)], with 
a scattering matrix M = M(ko) given by M = A^BA. 
Such relations can be extended to the case where the 
input waves g^ (x,0; k ) are not strictly plane waves, 
rather tilted broad beams with a narrow angular spec- 
trum G(An) = (27r) -1 / 2 j dxG(x) exp(— inx) and near- 
field distribution G(x). In this case we may write 

g (l) {x,0;K O ) = —= I cIAkG(Ak) x 

x exp ( ikqx + iAnx + i — Ix ) . (A7) 
V a J 

Repeating the previous analysis to each of the plane 
waves entering in the integral on the right hand side of 
Eq.(A7) and using the superposition principle yields 



g {l \x, z b ;k ) 



p=0,±l,±2, 



dAn—^=G{An)Mi ,<„ + As) exp + ,i« + 

V27T V a 



(A8) 



If the Fourier coefficients 9 nt i(n) of Bloch states - and spectral extension of G(An), we may take M.i :P {kq + 
hence the transfer matrix M(n)- vary slowly over the An) ~ M.i, p (kq) out of the integral in Eq.(A8). The re- 
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maining terms left under the integral then yields precisely which is extends the scattering matrix formalism to the 
g( p '(x, 0; Kq) [see Eq.(A7)]. We thus finally obtain case of broad beam excitation. 

p=0.±l,±2,... 

(A9) 
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